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ABSTRACT 

The Angstrom Project is undertaking an optical survey of stellar microlensing events across the bulge region 
of the Andromeda Galaxy (M31) using a distributed network of two-meter class telescopes. The Angstrom 
Project Alert System (APAS) has been developed to identify in real time candidate microlensing and transient 
events using data from the Liverpool and Faulkes North robotic telescopes. This is the first time that real-time 
microlensing discovery has been attempted outside of the Milky Way and its satellite galaxies. The APAS is 
designed to enable follow-up studies of M3 1 microlensing systems, including searches for gas giant planets in 
M3 1 . Here we describe the APAS and we present a few example light curves obtained during its commissioning 
phase which clearly demonstrate its real-time capability to identify microlensing candidates as well as other 
transient sources. 

Subject headings: gravitational lensing - galaxies: individual (M31) - techniques: photometric 



1. INTRODUCTION 

Gravitational microlensing has proved itself as a power- 
ful tool for planetary, stella r and Galactic studi es dPaczvhskil 
Il996t IBeaulieu et all 120061: iGould et al.ll2006l) . It has also 
been used to look for dark matter in the Andromeda Galaxy 
(M31) where the stellar crowding is much more severe 
dcalcni Novati et alj|2003t iRiffeser et alj|2003t Ide Jong et all 



12004 . Many of the exciting discoveries made by microlens- 
ing within our own Galaxy are possible only because mi- 
crolens ing is alerted in real time by survey teams such as 
OGLE (ISumi et al.l2006h and MOA (ISumi et al.l2003l allow- 
ing follow-up networks such MicroFUN ([Gould et al.ll2006l) 
and PLANET/RoboNet (IBeaulieu et al.ll2006l) to target can- 
didate events with high time resolution. The requirement of 
both rapid and robust data processing has until now prevented 
the development of an alert system for M31 microlensing. 

The Andromeda Galaxy Stellar Robotic Microlensing 
(Angstrom) Project aims to use microlensing to probe the 
geometry and ste llar mass function of the M31 bulge 
dKerins et al1 l2006). Since many of these events are expected 
to have durations of just a few days, the survey employs a 
distributed network of telescopes to enable round-the-clock 
monitoring of the M31 bulge. Two of the telescopes within 
the network are fully robotic, allowing data to be collected and 
processed rapidly. The Angstrom Project is taking advantage 
of this to develop a real-time data processing pipeline as part 
of a microlensing and transients alert system. The Angstrom 
Project Alert System (APAS) is designed to enable follow-up 
studies of candidate events in M3 1 in a similar vein to those 
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conducted within our own Galaxy. Chung et al.1 d2006l) has 
highlighted how such an alert system could enable microlens- 
ing discovery of gas giant planets in M3 1 itself. Other bene- 
fits of such a system would be a greater ability to identify and 
characterize binary microlensing events, wh ich are estimated 
to com prise around 15% of all M31 events dBaltz & Gondolol 
1200 lb . and the possibility of measuring finite source effects 
which can be used to constrain lens and source parameters. 

In this Letter we describe the design and implementation of 
the APAS and we show some early results from its commis- 
sioning phase during the 2006/7 Angstrom observing season. 
The results illustrate the potential of the system for real-time 
microlensing and transient discovery. 

2. ANGSTROM PROJECT OBSERVATIONS 

The Angstrom Project 8 employs a network of 2m-class tele- 
scopes. These include the robotic 2m Liverpool Telescope 
(LT) on La Palma, the robotic 2m Faulkes Telescope North 
(FTN) in Hawaii, the 1.8m telescope at Bohyunsan Observa- 
tory in Korea, the 2.4m telescope at MDM in Arizona and, re- 
cently, the 1.5m at Maidanak Observatory in Uzbekistan. The 
telescope network provides excellent round-the-clock cover- 
age of the M3 1 bulge during the observing season, which lasts 
from August through to February. Angstrom aims to obtain 
upwards of three observations per 24-hour period in order 
to be able to detect and characterize short-lived microlens- 
ing events t ypically induced by low-mass M dwarf stars in the 
M31 bulge (iKerins et alJl2006h . 

The APAS primarily utilizes incoming data from the two 
robotic telescopes, the LT and FTN. These telescopes operate 
autonomously, executing prog rams that ar e queued by an in- 
telligent scheduler program (lFraserll2006h . Images obtained 
by these telescopes are automatically pre-processed and are 
made available via the web typically within fifteen minutes of 
observation. Both telescopes employ identical optical cam- 
eras with a 4.6 arcminute field of view and 2kx2k CCD ar- 
rays. Angstrom observations are conducted primarily in the 
Sloan-like i band on both telescopes, typically with a se- 
quence of exposures totalling around 30 mins. Each expo- 
sure within the sequence is short (no more than 200 sees) in 

8 http://www.astro.livjm.ac.uk/angstrom/ 
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order to avoid saturating the core of the M31 bulge. The in- 
dividual exposures are then stacked to form a single 30 min 
ex posure image for eac h epoch. More details can be found 
in iDarnlev et al.l d2007l) . hereafter D07, which describes the 
main processing pipeline for the whole telescope network. 

In addition to LT and FTN data the APAS uses archive 
data from th e POINT- AGAPE dark m atter microlensing sur- 
vey of M31 dCalchi Novati et al.l2003l) . POINT- AGAPE used 
the wide-field camera of the 2.5m Isaac Newton Telescope on 
La Palma to survey a 0.6 deg 2 area of the M3 1 disk and bulge. 
Around 65% of the LT/FTN reference field overlaps with the 
POINT- AGAPE data. For the APAS we have reprocessed all 
Sloan /-band data from the POINT-AGAPE survey obtained 
between 1999 and 2001 . We use it to provide an extended data 
baseline for the subset of light curves lying in both Angstrom 
and POINT- AG APE survey regions. 

3. THE ANGSTROM PROJECT ALERT SYSTEM (APAS) 

M3 1 presents a challenging target for real-time microlens- 
ing discovery. At a distance of 780 kpc it is two orders of 
magnitude more distant than typical microlensing sources to- 
wards the Galactic bulge where the majority of Galactic mi- 
crolensing alerts occur. This means that larger telescopes and 
longer exposures are required to obtain sufficiently sensitive 
photometry. It also means that we have to deal with stellar 
fields which are far more crowded and where the typical mi- 
crolensing source is completely unresolved at baseline (see 
Figure[T|i. This is the so-called pixel-lensing regime. 

The challenge of developing an alert system for an M31 
pixel-lensing survey differs in some key respects from that of 
a Galactic microlensing alert system. One downside is that 
towards the core of the M3 1 bulge we reach the variable star 
crowding limit. Reliable photometry for individual variable 
obje cts becomes difficult even when using difference imag- 
ing dAlard & Luptonfl998l) . This means that very often more 
than one variable source contributes to the light curve at a 
given position, as illustrated in the variable source map shown 
in panel (c) of Figure Q] However, variable signals in the 
M3 1 bulge must necessarily involve a significant source flux 
change in order to be detectable above the very bright M3 1 
bulge surface brightness. Therefore there are fewer classes of 
variable stars that are sufficiently bright and variable to pro- 
vide potential false positive signals. 

With these potential difficulties in mind, we have adopted 
a pragmatic approach in developing the APAS. We do not re- 
quire it to autonomously and robustly flag up high quality mi- 
crolensing alerts. Rather the APAS presents a modest list of 
light curves that exhibit transient or microlensing-like behav- 
ior. The list is then reviewed by a human observer who makes 
the final decision on whether to issue an alert for any of the 
signals. The role of the APAS is to present to the observer 
via a web-based interface a compact list (typically no more 
than a few dozen light curves) of the currently most interest- 
ing transient signals from a continuously updated database of 
around 20000 light curves. The observer is then in a position 
to quickly review ongoing signals of potential interest. The 
APAS also records the positions of all variable signals in the 
immediate neighborhood of each light curve in the list so that 
the observer can check and quantify any light curve contami- 
nation. 

The APAS comprises two stages: a robust real-time data re- 
duction pipeline and an event classification scheme. We now 
describe both of these stages in more detail. 



3.1. Real-time data reduction 

The Angstrom Data Analysis Pipeline (ADAP) is respon- 
sible for reducing all Angstrom data from both robotic and 
non-robotic telescopes, and is described in detail in D07. The 
nominal mode of operation of the ADAP is to process the 
full dataset offline at the end of each observing season. How- 
ever we also employ the ADAP to process a lkx lk re-binned 
copy of the 2kx2k robotic telescope images that can be pro- 
cessed in real time. The pixel size of the re-binned data is still 
only 0."26 and so adequately samples the typical seeing disk 
(1" - 1 ."5). The ADAP offline processing is described fully in 
D07, whereas here we concentrate on the real-time processing 
of the re-binned data which forms part of the APAS. Many of 
these stages are similar or identical to the ADAP offline pro- 
cessing so we keep our description brief and refer the reader 
to D07 for further details. 

Firstly, an initial quality control score is made of each frame 
based upon the shape and size of the point-spread function 
(PSF). Images which fail this check are dropped from the sub- 
sequent analysis as their inclusion would result in unaccept- 
able difference image quality. For the surviving images the 
APAS performs defect masking, cosmic ray rejection and im- 
age alignment. Then the APAS de-fringes each frame using 
an iterative method described in D07. Images taken within 
the same run (i.e sequentially at one epoch) are then PSF 
matched and stacked using the AngstromlSIS difference im- 
age engine, which is a modified version of a s ubset of routines 
from the ISIS image subtraction package dAlard & Luptonl 
119981) . AngstromlSIS is described in detail in D07. The im- 
age stack selected to be the reference is then convolved and 
subtracted from the other image stacks to create a sequence of 
difference images. These are then used to produce one likeli- 
hood map of variable sources for each epoch, using an indi - 
cator of statistical significance similar to that of Cash ( ' 1979b . 
Sources from within each likelihood map are then identified 
and matched and an updated master list of variable sources is 
compiled. PSF-fitting photometry is performed on the differ- 
ence images and the resulting photometry is stored within a 
MySQL database. The time between observation and adding 
new photometric points to the database is about two hours. 

This pipeline is applied both in real time to the LT and 
FTN data and offline on regions of legacy POINT-AGAPE 
data which overlap with the Angstrom field. The LT and FTN 
are essentially identical telescopes and the observations were 
taken using matched filter systems, so a simple linear scaling 
is sufficient to deduce relative light curve scalings and off- 
sets for these telescopes. The POINT-AGAPE /-band data is 
matched to LT/FTN /-band photometry by fitting a sample of 
periodic variables present in both the Angstrom and POINT- 
AGAPE fields. 

3.2. Candidate classification 

APAS event classification is initiated each time a new 
robotic observation is added to the light curve database. The 
APAS employs a number of selection criteria to classify each 
light curve into one of seven different categories. 

Each object's light curve is separated by telescope (LT, FTN 
or POINT-AGAPE), with each telescope being initially pro- 
cessed separately. The APAS deems an object to be valid if 
it contains at least seven LT light curve points from the first 
(2004/5) and second (2005/6) Angstrom observing seasons. It 
must also contain at least seven FTN light curve points from 
the second Angstrom season. If an object does not meet these 
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FIG . 1 . — (a) LT ('-band image of a two arcminute region of the Angstrom M3 1 bulge field, (b) Corresponding difference image obtained after convolving and 
subtracting a reference image, showing variable objects as black and white spots. The residual of the imperfectly subtracted core of the M31 bulge is visible in 
the top right-hand corner, (c) The corresponding significance map of variable sources visible above the background. Variable source blending is evident across 
the map. Blackened regions denote masked areas, (d) Stars within a two arcminute region of an OGLE-III Milky Way bulge field shown for comparison with the 
crowding levels seen in the variable source map in (c). Whilst all objects in (c) have variable flux only a handful of the stars visible in (d) are variables. OGLE 
image is from the OGLE Early Warning System lUdalski 2003) for alert OGLE-2005-BLG-172. 



criteria it is classed as Void and will never be processed again 
by the APAS. If a light curve does not contain at least seven 
points in the current (2006/7) observing season then the object 
is marked as Invalid. Such objects are reviewed again when 
new observations are added to the database. 

For remaining objects a baseline is estimated for each tele- 
scope by determining the median light curve flux. The APAS 
then looks for indications of transient behavior. Any point that 
is at least 5er above the baseline is considered "significant". 
For each object the APAS then investigates any clusterings 
of significant points ("peaks") for each telescope. A peak is 
defined as a clustering of at least five significant points, all ly- 
ing within the same observing season. Within the peak some 
leeway is permitted - up to three consecutive non-significant 
points are admitted within the peak provided they are fol- 
lowed by a significant point within the same season and pro- 
vided their presence coincides with worsening seeing. 

Objects containing no peaks are classed as Flat. These 
objects will be re-analyzed and possibly re-classified when 
new observations are added to the database. If an object con- 
tains no peaks in the current season but contains at least one 
peak within a previous season then it is classed as Variable, 
undergoing re-analysis when additional observations are ob- 
tained. Objects containing peaks in the current season are 
checked for the presence of peaks in previous seasons. The 
APAS must allow for the fact that the crowding of vari- 
ables may superpose variability on the baseline of transient 
and microlensing events (for example candidate PA-99-N1 
iPaulin-Henriksson et al.l l2003). Accordingly, it does not nec- 
essarily classify objects as variable if they show variability 
on previous seasons as well as in the current one. Instead, if 
peaks are seen in previous seasons then the object is classified 
as variable unless a peak in the current season is at least 150 
reference image counts and is at least 50% greater than the 
highest amplitude peak in any previous season. Other than 
these the only remaining objects are those which show only 
peaks during the current season. 

The APAS now considers all surviving objects to be tran- 
sients and flags them as Followed. The FTN and POINT- 
AGAPE (where available) light curves for these objects are 
now scaled to the reference image of the LT data and a re- 
duced Paczyhski profile is then fitted with a flux F(t ) = B + 
AF I yl + 12[(f-fo)/fFWHivi] 2 , where t and fn are the epochs 
of observation and maximum brightness, B is the baseline 
flux, AF = F(to)-B is the maximum flux deviation and ?fwhm 
is the duration of the full-width at half-maximum. The values 



Table 1 

a snapshot statistical summary from the apas. 



Object 


Number of Percentage 


status 


objects 


of total 


Void 


6,311 


34.0% 


Invalid 


54 


0.3% 


Flat 


9,931 


53.6% 


Variable 2,197 


11.8% 


Follow 


12 


0.1% 


Old 





0.0% 


Alert 


38 


0.2% 



of B and AF are telescope and bandpass dependent. We re- 
quire that the best fit to does not occur before the first obser- 
vation of the current season and that it does not occur after the 
anticipated start of the next season. The FWHM of the event 
must be in the range of 1 day < ?fwhm < 365 days. Finally 
we require that the amplitude of the event be positive. Ob- 
jects passing these criteria are classed as Alerts. Objects that 
fail any of these criteria remain classed as Followed, unless 
they were previously classified as an alert in which case they 
are re-classed as Old. 

The APAS utilizes a web-based interface to report results. 
Whilst allowing a user to examine the light curves of all 
alerted and followed objects it also displays light curves of all 
objects within a three arcsecond vicinity of each alert. This is 
helpful in deciding whether an object is a true alert but may 
have a baseline contaminated by a nearby variable source. 

3.3. Results from the commissioning phase 

The Angstrom Project is currently in its third observing sea- 
son and this is the first season in which we have implemented 
the APAS. Table Q] contains a snapshot summary of statistics 
from a recent run of the APAS. Around 12% of the light curve 
catalogue is classified as variable or transient. The bulk of the 
rest of the light curves are either "Flat" because they do not 
show significant light curve variation, or "Void" because they 
do not comprise enough data points in previous observing sea- 
sons. Reassuringly, the number of variables flagged up by 
the APAS is in line with the variable number density derived 
from an analysis of thre e years of data from the wide area 
POINT-AGAPE survey (lAn et al.ll2004b . This indicates that, 
whilst probably not optimal, the APAS nonetheless competes 
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FIG. 2. — Three examples of alerts selected by the APAS. LT photometry is shown in black and FTN in red. (a) A short time-scale (fFWHM — 1 day) high signal- 
to-noise ratio microlensing candidate, (b) A lower signal-to-noise ratio microlensing candidate at the borderline of identification by APAS . Both micro lensing 
candidates exhibit periodic variability in their baselines due to nearby variable sources, (c) Classical Nova 2006 8" first announced by Burwitz et al. i200g). This 
light curve is easily distinguishable from microlensing by its asymmetric profile. Right-hand panels show zooms of the events and the solid line is the best-fitting 
reduced Paczyriski curve. 

"Rochester naming convention http://www.rochesterastronomy.org/snimages/index.html 



in terms of sensitivity with previous non real-time pipelines. 
We can expect that the offline ADAP analysis of the un-binned 
data will be more sensitive still. APAS also succeeds in pre- 
senting an easily manageable number of alerted and followed 
light curves for a human observer to review. 

Three light curves from the subset of 38 listed alerts are 
shown in Figure [2] Clearly, these examples represent credi- 
ble microlensing and nova candidates. The APAS therefore 
demonstrates the capability to provide a rapid, robust and 
credible extragalactic alert system. The first candidate in par- 
ticular appears to be a very short high signal-to-noise ratio 
microlensing event. Such events tend to be sensitive to finite 
source effec ts and this event i s currently the subject of a sep- 
arate study (Du ke et al. II2007I) . 

4. DISCUSSION 

We have successfully implemented the Angstrom Project 
Alert System (APAS), which performs real-time reduction 
and identification of extragalactic transient sources and mi- 
crolensing events. This is the first system of its type to pro- 
cess extragalactic observations. The APAS has been com- 
missioned during the third season of Angstrom observations 
and is already displaying the capability to identify credible 



microlensing and nova candidates. During the remainder of 
the current 2006/7 observing season and in the following ob- 
serving seasons we aim to use the APAS to flag up candidate 
sources for higher cadence follow-up in order to obtain de- 
tailed light curves of novae outbursts, and to identify exotic 
microlensing systems. 
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